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Abstract. The standard assumption for the equilibrium microcanonical state in 
^ I quantum mechanics, that the system must be in one of the energy eigenstates, 

is weakened so as to allow superpositions of states. The weakened form of the 
, microcanonical postulate thus asserts that all quantum states giving rise to the same 

energy expectation value must be realised with equal probability. The consequences 
that follow from this assertion are investigated. In particular, a closed-form expression 
^ I for the density of states associated with any system having a nondegenerate energy 

' spectrum is obtained. The result is applied to a variety of examples, for which the 

I . behaviour of the state density, as well as the relation between energy and temperature, 

' are determined. Numerical studies indicate that the density of states converges to a 

■ distribution when the number of energy levels approaches infinity. (1 February 2008) 

o 
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Q^, 1. Introduction 

' is sometimes argued that the quantum-mechanical description of the microcanonical 

' distribution is simpler than the corresponding classical counterpart 113]. This is 
because the standard quantum microcanonical postulate asserts that if the energy of 
the system lies in the range E to E + AE, then all of the energy eigenstates with energy 
En G [E,E + AE] are realised with equal probability Pj. While such a postulate does 
indeed provide an elementary statistical description of the quantum system in terms of 
a uniform probability distribution, it is not clear how this postulate ties in with other 
fundamental ideas in quantum mechanics. 

For example, in quantum mechanics the state of the system can be represented by a 
general superposition of energy eigenstates, whereas according to the standard quantum 
microcanonical postulate, superpositions of energy eigenstates with distinct eigenvalues 
are excluded. Also, if an energy eigenvalue that lies in the interval [E,E + AE] is 
nondegenerate, then the standard postulate implies that the system can only be in 
a single eigenstate corresponding to that eigenvalue. These observations suggest that 
the standard postulate is perhaps too stringent to give rise to a satisfactory statistical 
description of a quantum system in isolation. 
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The purpose of the present paper is to study the consequences of a relaxation of 
the standard quantum microcanonical postulate. Specifically, we consider the following 
generalisation of the microcanonical postulate: namely, that every quantum state 
possessing the same energy expectation value must be realised with an equal probability. 
According to this 'weaker' postulate various superpositions of energy eigenstates with 
distinct energies are no longer excluded; likewise, a microcanonical state will never 
correspond to a single eigenstate (except in the case of the smallest and the largest 
energy levels if these are nondegenerate) . 

The idea that we investigate here is similar to the classical case, although there 
is a subtle difference. Classically, the uncertainty in energy is fully characterised 
by a statistical distribution over the phase space, and for a classical microcanonical 
distribution having support on a level surface of the Hamiltonian the energy variance 
vanishes. Quantum mechanically, the contribution to the energy variance from the 
statistical distribution over the phase space also vanishes. However, if the specified 
energy level is not the largest or smallest energy eigenvalue, then there is an additional 
contribution to the energy variance that arises from the intrinsic quantum uncertainty. 
Therefore, according to our quantum microcanonical postulate the energy variance of a 
system described by a quantum microcanonical ensemble need not vanish, even though 
the system is in isolation. 

The paper is organised as follows. In ^ we review briefly the formulation of 
standard quantum theory in terms of the geometry of the space of pure states, as 
described, e.g., in [H El El El El CHI references cited therein. When the trajectory of 
a wave function is projected from Hilbert space Ti. to the space V of pure states (rays 
through the origin of H), the Schrodinger equation on H reduces to Hamilton's equation 
on V. More precisely, the space of pure states V has a natural symplectic structure; 
and the Schrodinger trajectories, when projected onto V, are the integral curves of the 
Hamiltonian vector field obtained by taking the symplectic gradient of the function 
on V defined by the expectation of the Hamiltonian operator. As a consequence, the 
formulation of quantum mechanics on V provides a natural environment in which one 
can study important issues arising in the context of Hamiltonian mechanics, such as 
ergodicity conditions or the construction of equilibrium ensembles. For example, the 
dynamical approach to microcanonical equilibrium introduced by Rugh ^2] for classical 
systems can be seen to apply in the quantum regime. 

In ^ we introduce the microcanonical density of states that follows from our 
weakened microcanonical postulate, together with the associated microcanonical density 
matrix. In 21 "we derive a general integral representation for the density of states, 
expressed in terms of the energy eigenvalues. We then perform the integration explicitly 
in ^in the case for which the energy spectrum is nondegenerate. As an illustration, the 
properties of a system having an equally spaced energy spectrum are studied in detail 
in ^ In this example we study the relation between the energy and the temperature 
in some detail, which we plot for a number of different situations. In particular, a 
simple procedure for rescaling the energy shows that as the dimensionality of the Hilbert 
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space increases the system becomes more and more likely to take the intermediate 
energy value ^(-Emax + -Emin), where -Emax and -Emin are, respectively, the largest and the 
smallest energy eigenvalues. To study the convergence of the distribution numerically, 
we compute in 33 the Hellinger distance between the density of states associated with 
an n-level system and an (n + l)-level system for n = 2, 3, . . .. The result shows that 
the logarithmic plot of the relative distances of the distributions against the Hilbert 
space dimensionality lies on a straight line with gradient —2. Finally in ^ we study 
the properties of systems having other nondegenerate spectra. In particular, we show 
numerically that there is an approximate symmetry relation that holds between a system 
having the spectrum En ~ n'^ and a system having the spectrum En ~ n^/'^, where k is 
a constant. 

2. Quantum phase space 

In order to investigate properties of a closed, isolated quantum system in equilibrium it 
will be useful first to recall the Hamiltonian formulation of standard quantum mechanics. 
This has the advantage of allowing us to apply concepts arising in the corresponding 
classical theory of equilibria as outlined, for example, in Ref. 0. 

Physical states in quantum mechanics are represented by elements of a complex 
Hilbert space 7i, which we assume to be (n + l)-dimensional. Let us denote by 
a typical element of TC, so that the index a runs over the range a = 0, 1, . . . ,n. The 
Hamiltonian, which acts on elements of Ti, can thus be denoted Hj^, and the expectation 
value of Hp in the state Z'^ is given by the expression 

where Za is the complex conjugate of Z". The Hilbert space Ti. carries an essentially 
irrelevant complex degree of freedom given by the overall scale of the state vector. This 
follows from the fact that the expectation value of a physical observable is invariant 
under the complex scale transformation Z°' —>■ AZ", where A G C — {0}. It is useful 
in some applications to eliminate this extra degree of freedom by considering the space 
of equivalence classes under the relation XZ°' ~ Z°' for A G C — {0}. This is the space 
of rays through the origin of Ti, otherwise known as the projective Hilbert space V 
of complex dimension n. It is well known that quantum theory, when formulated on 
the projective space V, admits a representation in terms of the standard mathematical 
structure of Hamiltonian mechanics. This can be seen as follows. 

We find it convenient for our purposes to regard the projective space P as a real 
manifold F of dimension 2n, letting (a = 1, 2, . . . , 2n) denote a typical point in F. 
Therefore, each point x*^ G F represents a ray in the Hilbert space Ti. In this way we 
can regard the expectation as determining a real-valued function H{x) on F. 

The space of pure states, when regarded as the real even- dimensional space F, is 
endowed with a symplectic structure, given by a nondegenerate, antisymmetric two- 
form Uab- The dynamical laws governing the trajectories of quantum states, given by 
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the Schrodinger equation on Ti, can then be represented on T in Hamiltonian form as 
follows: 



In other words, the space F is a symplectic manifold upon which the evolution of a 
quantum state is governed by Hamilton's equations, which in the language of symplectic 
geometry take the form 0- Therefore, we can regard F as the quantum analogue of 
a classical phase space. As a consequence we can also formulate our investigation of 
equilibrium states on F. 

3. The microcanonical ensemble 

We begin this section by considering the foliation of the quantum phase space F by level 
surfaces of the Hamiltonian function H{x). This is given by a family of hypersurfaces 
{Se}, E G [i^min, -Emax], determined by level values H{x) = E of the Hamiltonian 
function. The structure of the typical energy surface in quantum mechanics is quite 
intricate, even for a system described by low- dimensional Hilbert space. In particular, 
as E varies in the given range [Emm, -f'max] both the dimensionality and the topology 
of the associated energy surfaces can change. An example can be found in Ref. [HI, in 
which the structures of the energy surfaces for a three-level system are investigated in 
detail. 

Now given this foliation, the 'number' of quantum mechanical microscopic 
configurations (pure states) with expected energy in the small range E and E + AE is 
fl{E)AE, where the state density Q{E) for energy E is given by an expression of the 
form 



where ebc-.-d denotes the totally skew tensor with n indices, and 

gab ig 

the inverse of the 

natural Riemannian metric Qab on F. The metric Qab is compatible with the symplectic 
structure Uab in the sense that VaUJbc = 0, where Va is the unique torsion- free covariant 
derivative operator on F satisfying 'W aQbc = 0. It is a remarkable feature of the quantum 
phase space F that it has both a natural Riemannian metric and a compatible symplectic 
structure. These elements can be regarded as part of the natural geometry of any 
quantum system. 

In the case of an isolated quantum mechanical system with energy in the small 
range E to E + AE, we can adopt the notion of the microcanonical ensemble in classical 
statistical mechanics (cf. [131 )j identify the entropy of the system by use of the 
Boltzmann relation 




(2) 




(3) 



Here the natural vector- valued {2n — l)-form da"' on F is defined by 
da'^ = (?"V..ddx'=---dx^ 



(4) 



S{E) = kln{n{E)AE). 



(5) 
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Here we are implicitly assuming what might be called the quantum microcanonical 
postulate, which asserts that for an isolated system in equilibrium all states on a given 
energy surface in the quantum phase space are equally probable. As a consequence, the 
temperature T of such a system is determined by the relation 

where f3 = 1/kT and k is Boltzmann's constant. Thus for an isolated quantum 
system with expected energy E, we assume that the equilibrium configuration is given 
by a uniform distribution on the energy surface Se, with entropy S{E) and inverse 
temperature l3{E), as given above. The corresponding probability density on F, which 
we call the microcanonical F-distribution, is 

f^Eix) = ^^6{H{x)-E), (7) 

where 

n{E) = j 5{H{x) - E)dV. (8) 

Here dV is the volume element on F. 

It is a straightforward exercise to show that starting from the definition (jHJ of the 
density of states Vt{E) we can deduce the integral formula This can be seen as 
follows (cf. dH). 

First we note that at each point x G F such that Va-ff(x) 7^ the volume element 
on F can be written as a product 

<W = dNda, (9) 

where the (n — l)-form da defined by 

, V''Heab-cdx''---dx^ 

da = , ^ = (10) 

is the volume element on the energy surface passing through x, and 

dN= ^-'^^^ (11) 

On the other hand, as a consequence of the relation 

dH = VaHdx" (12) 
we observe that 

dN = ^ ^ . (13) 



Substituting this expression into (jH)) we obtain 
n{E) = j 6{H{x) - E)dV 



j 6{H{x) - E) 



dadH 
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V'Heab-cdx' ---dxl ^^^^ 



which agrees with expression (jS)). 

A general measurable function F{x) on T represents a nonlinear observable in 
the sense of Kibble ^U] and Weinberg The usual 'linear' observable of standard 
quantum mechanics then corresponds to the situation for which F{x) can be represented 
as the expectation of a Hermitian operator. In either case, for each value of x we interpret 
F{x) as the conditional expectation {F)x of the observable F in the pure state x. The 
unconditional expectation of F in the microcanonical F-ensemble is then given by 



{F)e = J^Fix)fiEix)dV. (15) 

In the case of a linear observable we have 

Fix) = Fi*n^(x), (16) 



where 



n"(x) = — (17) 

is the projection operator onto the state vector Z"{x) corresponding to the pure state 
X G r. Then the unconditional expectation in the state ^e{x) is 

{F)e = F^f^-piE), (18) 

where the quantum microcanonical density matrix fi'^{E), parameterised by E, is defined 
by 



/i^(E) = / n^(x)/x^(x)dK (19) 



r 



Providing we only consider linear observables, i.e. as in standard quantum mechanics, 
the state of the system is then fully characterised by the density matrix fi'^{E). 

Now suppose W{E) denotes the total phase space volume for states such that 
H(x) < E. Then the density matrix fJ^'^iE) can be calculated explicitly by use of a 
'variation-of-parameters' formula given by 

..^)^-(^)-^. 

This representation can be verified as follows. From the definition of W{E) we can write 

WiE) = r [s (i/f n^(x) - u) dVdu, (21) 
J-oo Jr 

and hence 

-E r 

Il'^ix)6' {H^W^ix) - u) dVdu 



dW{E) 



dm 



oo Jr 
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n^(x)5(i7fn^(a;)-E) dV 



= -n{E) J n°(x)/is(x)dr 

= -niE)^^^{E). (22) 
On the other hand, clearly dW{E)/dE = fl{E), and thus we obtain (PU)). 

4. Calculating the density of states 

As defined in equation (jH)) above, the density of states is given by the volume integral 
over r of a delta-function having a support on the energy surface Se- Our objective now 
is to perform the relevant integration explicitly for a generic Hamiltonian, and obtain a 
representation for Q{E) in terms of the energy eigenvalues. 

We find it convenient to pursue the calculation by lifting the integration from the 
phase space F to the Hilbert space Ti, imposing the constraint that the norm of the 
Hilbert space vector is unity. Therefore, we write the expression (jSI) in the following 
form: 

Vt{E) = - / - 1) 5 E d"+iZ d"+iZ. (23) 

The additional factor of vr arises here from the superfiuous phase integration in 
By use twice of the standard integral representation 



1 



oo 



S{x) = / e-'^MA (24) 



we thus deduce that 



7r ./_^ 2-K 2ti 



X [ exp [(-i (AZ^Z" + vH^Z^Z^) 1 d^+^Z d"+^Z. (25) 

We now observe that we can diagonalise the Hamiltonian by unitary transformation 
without affecting any of the terms in ()25p on account of the fact that every 'kef vector 
Z" is coupled to a corresponding 'bra' vector Z^. Therefore, the density of states can 
be written in the form 

n{E) = ^ r l^e^(A+^^) f d-^^z d-^'Z exp ( -i V(A + uE,)Z,z] , (26) 

where {-Ei}/=o.i....,n are the energy eigenstates. 

This is of course a formal expression; the integration can nevertheless be carried 
out if we regard ()26|) as the limit of a similar integral in which A and u are displaced into 
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the complex along the negative imaginary axis. The integration over C""*"^ then reduces 
to a {2n — 2)-dimensional Gaussian integral, which is readily performed to yield 

This is the desired integral representation for the density of states, expressed in terms 
of the energy eigenvalues. 

5. Density of states for a nondegenerate energy spectrum 

We proceed further by evaluating the integration in (j27p in the case where the 
Hamiltonian has no degenerate eigenvalues. Let us consider the integration in the A 
variable first. We observe that there are n + 1 first order poles on the real A- axis. An 
application of Cauchy's theorem thus gives us 

from which it follows that 

Q{E)=n^^y / ^iJ^ TT —1--. (29) 

We now recognise the //-integration formally as the ra-fold repeated integral of the 6- 
function, defined by the truncated polynomial 

*'-">w-{"^.- til <^°' 

As a consequence, the density of states associated with a quantum system having a 
nondegenerate energy spectrum can be seen to be given by an expression of the form 
given by 

n n ^ 

n{E) = {-i)vJ2^'~-\E,-E) n Y-^- ^2^) 

k=0 l=Q,j^k ' ^ 

In addition to the function Q{E) we find it useful for some purposes to introduce 
the related normalised density of states fi{E) defined by 

m ^ . (32) 

The function /i(-E') thus satisfies the property that 

/oo 
fi{E)dE = 1. (33) 
'OO 

To calculate the normalisation factor we need to determine the total volume Vr of the 
quantum phase space. This is given by 

Vr = - [ 5{Z^Zo, - l)d"+^Zd"+^Z, (34) 

TT 7r"+i 
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where again the factor of vr in the denominator refers to the removal of the superfluous 
overall phase of the wave function. Thus, making use of the integral representation of 



the delta-function, we have 



Vr = - [ d^+'Zd^+'Z r —e'^i^^^^-'^dX 

1 rdX^^_,^r giAZ'^Z.^n+l^^n+l^ 



-iX 



i / aA :i / ITT 



— e 

— , 35 
nl 



where in the last step we have used the identity 
1 e-'^ , , i" 



dA = -r. (36) 

This is in agreement, e.g., with the result obtained in Ref. |7j by use of different methods. 
In particular, we note that in the case n = 1 (a two-dimensional Hilbert space) the space 
of pure states T is isomorphic to the geometry of a round two-sphere with radius one-half. 

6. Equally spaced energy spectrum 

We now specialise the analysis further by considering a quantum system whose 
energy eigenvalues are equally spaced. There are many elementary systems, such as 
noninteracting particles with spin, having a spectrum of this kind. To fix the units of 
energy we set 

Ek = hujk, (A; = 0,l,...,n). (37) 
Substituting this in (pTH) we obtain 

^ tll^^f^^yn-i y ff —1—. (38) 

^ ' k>[E\ ^ ^ l=0,^k ^ ' 

This expression can be simplified further by writing 

1 _ 1 
nuj{l - k) ~ {huj)''{-l)^k\{n- k)\ 



l=0,^k 



and hence the microcanonical density of states associated with a system having equally 
spaced energy spectrum is given by 



\n-l 



^ ^ k>[E] ^ ' 

Finally, dividing this expression by the normalisation factor (j35|) we obtain the following 
formula for the normalised density of states: 

k->{E\ ^ ' 
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It will be useful now to indicate what this function looks like in various situations. In 
Figure we plot fi{E) for two, three, and four dimensional Hilbert spaces. 

KE) f'S) MB 




-1 -0.5 0.5 1 1.5 2 -1 -0.5 0.5 1 1.5 2 -1 -0.5 0.5 1 1.5 2 

Figure 1. Density of states as a function of energy E for systems witlr two, three, 
and four nondegenerate energy eigenstates. Tlie functions are evaluated piecewise. 
For example, in the case of the four-dimensional system with En/hw = 0, 1,2,3, the 
normalised density of states ^l{E) is zero for E <Q and E > i, and is given by three 
distinct quadratic functions in the intervals (0, 1], (1,2], and (2,3]. In general, for an 
n-dimensional system, ij.{E) is given by a combination of polynomials of degree n — 1, 
and is at least n — 3 times differentiable for all values oi E. In all cases the area under 
the function integrates to unity. 




Figure 2. Rescaled normalised density of states /i(-E) as a function of energy E. As we 
increase the density of energy levels in the interval [0, 1], the corresponding density of 
states becomes more sharply distributed around the centre E — ^. For these numerical 
plots we have set fno — \. The plots correspond to systems for which the number of 
energy levels ranges from 3 to 12, and the value of the rescaled energy varies over the 
unit interval. 



We would like now to study the behaviour of the density of states as we increase 
the number of energy levels. For this purpose we find it convenient to rescale the energy 
spectrum so that the range of energy is over a fixed interval [0, 1]. After the application 
of a suitable such rescaling, the density of states reduces to 

k>\E] ^ ' 
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In Figure 121 we plot this rescaled density of states for a variety of systems with different 
numbers of energy levels. 

Expression for the density of states also allows us to study the relation between 
energy E and temperature T, by use of the relation In Figure El we plot the system 
energy E{T) as a function of temperature by numerically inverting the relation 

T^{§y\ME). (43) 

As the temperature increases from zero, the corresponding energy increases 
monotonically, and asymptotically approaches the value = | (in the rescaled units). 
The function E(T) is multi- valued at T = 0, where the energy takes the values E = 0,1. 
In the region E G (|, 1] the density of states fi{E) is a decreasing function of E, and 
thus the corresponding 'temperature' is negative. To put the matter differently, this is 
the region that is not accessible in an equilibrium. 



E(T) 




0.1 0.2 0.3 0.4 0.5 

Figure 3. Energy E{T) as a function of temperature T. The linear curve with lowest 
energy corresponds to a four-level system, and we have plotted E(T) up to the ten-level 
system with an interval of 1. Although not shown in the plot, the second root for E{0) 
is given by E(0) = 1, where the function is multi- valued. As T is reduced from zero, 
E{T) reduces monotonically, and approaches E — ^ as T ^ — oo. This is the region 
inaccessible in equilibrium. 



7. Convergence in the infinite energy-level limit 

In statistical mechanics it is often presumed that the various different distributions 
used to describe equilibrium states (e.g., microcanonical, canonical, grandcanonical, 
and pressure-temperature distributions) should in some respects be equivalent in the 
infinite volume or thermodynamic limit. Now the system that we have studied in the 
previous section corresponds to a system of noninteracting quantum particles. As a 
consequence, we would not expect any nontrivial behaviour exhibited by the system 
in the thermodynamic limit. Nevertheless it would be of interest to study how the 
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density of states behaves in this hmit so that the result might be compared with the 
corresponding canonical formulation outlined in Ref. 0|. 

Now the results shown in Figure El indicate that as we increase the number of energy 
levels the density of states becomes more and more peaked at the intermediate energy 
E = ^. Therefore, in this section we consider the separation of a pair of normalised 
state densities associated with /c-level and {k + l)-level systems for a range of values for 
k, and study whether the separation decreases as we increase k. If so, then the result 
will indicate that the density function is converging, possibly to a delta-function centred 
at E = i. 

There are various standard measures that one can use to study the separation of 
a pair of density functions, such as the relative entropy or the Bhattacharyya distance. 
Here we shall consider a closely related measure given by 



(44) 



known as the Hellinger distance. For clarity we let fik{E) denote the density of states 
associated with a fc-level system. Note that the Hellinger distance is simply the L^-norm 

D(/i„,/i„)= ^/JiJ^ - ^/J^Je) (45) 

of the difference between the two square-root density functions. 



Log D(Aik,//k+i) 




Log k 



Figure 4. The Hellinger distance D{iik, tJ-k+i) between two 'consecutive' rescaled state 
densities is plotted for fc = 2, 3, . . . , 10. The plot is presented on a logarithmic scale. 
We find that the logarithms of the distances lie almost perfectly on a straight line with 
gradient of —2, indicating that the separation of the densities decays quadratically 
with an increasing number of energy levels. 

The result of the numerical analysis of ()44j) is shown in Figure 01 which indicates 
that the relative separation of fik{E) and ^k+i{E) is decreasing quadratically in k. 
Therefore, we conclude that the density of states is likely to be converging in the sense 
of L2. 
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8. Other spectral structures 

In the example considered in ^Hl and ^ above we analysed a linear spectral structure 
of the form oc k, and found that the density of states is centred at an intermediate 
energy value. In this section we investigate nondegenerate systems having other kinds 
of energy growths. We start with the case of a system whose energy grows quadratically 
so that El. = hwh"^. In this case we make use of the relation 



n 



n 



1 



l=0,^k l=0,^k 

and substitute this in formula (j^T|) to obtain 

2ni-ir A (-1)' {k^ 



{l + k){l-k) {n + k)\{n-k)\ 



(46) 



n 

E 

k>[E] 



)' 



[n + k)\{n - k)\ 



(47) 



As one might have expected, the functions given by ()47|) have similar characteristics 
to those of the previous example. For each value of n we obtain a curve fi{E) which 
is specified by n + 2 polynominals of degree n — 1. The function is continuous and its 
differentiability is of order n — 2. The first three examples are plotted in Figure El As 
the plots indicate, the microcanonical density functions in the current examples are no 
longer symmetric around the intermediate energy |(-Ejnin + -Emax), in contrast to the 
case of linear energy growth. 





Figure 5. The density of states Hk{E) for fc-level systems witli fc = 2, 3, 4, in tlie case 
of a system with a quadratically growing energy spectrum. For the plot we have set 
hw — 1. Unlike the example with a linear, equally-spaced spectrum, the density of 
states is not symmetric. 



As in the previous example, in the present case we can rescale the energy spectrum 
in such a way that we can directly compare the behaviour of the density of states as we 
increase the number of energy levels. The result is shown in Figure El We observe that 
the density of states becomes more peaked as we increase the number of energy levels. 
However, the location of the peak is no longer at = |, but rather closer to -E = |. 

As a consequence of the skewed form of the distribution, the range of energy for 
which the derivative of ^i{E) with respect to E remains positive is somewhat reduced. 
This implies that the range of energy associated with positive temperature is reduced 
from the previous example of linear energy spectrum. Some examples are illustrated in 
Figured where we plot the temperature dependence of the energy. As the temperature is 
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increased, the energy grows monotonically and reaches a value around i. The remaining 
values of the energy are associated with negative temperatures that are inaccessible in 
equilibrium. 



ME) 




Figure 6. The density of states ^i{E) associated with a system having the spectrum 
Ek = k'^/{N—l)^, where N denotes total number of energy levels. The value of k thus 
ranges over = 0, 1, . . . , — 1. Here we plot /u(-E') for A'' = 4, 5, ... , 11. Numerical 
studies show that the peaks of the distributions are located at approximately E = ^. 



We have examined systems having a linear energy growth and a quadratic energy 
growth. In the linear case the microcanonical distribution is symmetric around its 
centre, while in the quadratic case the peaks of the distributions have shifted to the 
left with smaller energies. This leads to the question of what happens to the density of 
states associated with systems having other spectral structures. 




0.5 



Figure 7. System energy E{T) as a function of temperature T in the case of a system 
having a quadratic energy spectrum. The plots correspond to a set of fc-lcvcl systems 
with k = 5, 6, . . . , 11. Although not shown in the plot, the energy is multi-valued at 
T = so that E{T) ^ as T ^ 0+ and E{T) ^ 1 as T ^ 0". The range of energy 
accessible in equilibrium is thus given hy < E < i. 
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In the case of a system with an infinite number of degrees of freedom (and thus an 
infinite number of energy levels), the energy spectrum cannot grow more rapidly than 
quadratically in the number of energy levels. However, for a finite system there is in 
principle no limitation on how fast the system energy can grow. Therefore, we would 
like to study the behaviour of the density of states associated with systems having the 
following two spectral structures: 

E„ = and E„ = n^/^ (48) 

for k = 1,2,3,.... Thus in one case the growth of the energy is enhanced as k is 
increased, while in the other case the growth of the energy is suppressed. Various 
densities of states associated with systems having spectra of the form PHj) are computed 
numerically and plotted in Figure |H1 The result shows that as the rate of growth is 
increased, the location of the peak of /i(-E) becomes smaller, whereas when the rate of 
energy growth is suppressed, the peak of yu(-E) increases. In particular, the peak of fi{E) 
for the system with energy En = n'' is located approximately around E ~ (A; + 1)~^, 
while the peak of n{E) for the system with energy En = n^^^ is located approximately 
around E ~ 1 — (A; + l)~^. Therefore, when /c <^ 1, the values of energy that are 
accessible in equilibrium becomes negligible when E^ = whereas li E^ = ut^ then 
virtually all values of G [0, 1] are accessible. 



ME) 




Figure 8. The density of states /i(i?) associated with systems having the spectra 
En — and En — , for fc = 1, 2, . . . , 10. In the former case the peak of n{E) shifts 
to the left, while in the latter case the peak shifts to the right. The cases £"„ = and 
En = nk are close to but not exactly symmetric to one another around E = ^. In all 
examples the number of energy levels is set to six. 

In conclusion, we note that the analysis we have pursued here is based on 
a generalised microcanonical postulate, which relaxes the somewhat more stringent 
assumptions made in the standard quantum microcanonical postulate as outlined in, 
e.g., Ref. jH|. To determine whether the assumption made in this paper refiects the 
actual equilibrium distribution of an isolated quantum system it will be desirable to 
examine the properties of interacting systems. In the case of a general interacting 
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system, some of the energy levels are typically highly degenerate, and thus we must 
return to the integral representation (jTTj) with a view to deriving an efficient way to 
carry out the integration in the cases when there are higher-order poles. This is an 
intriguing open problem that we hope to investigate elsewhere. 
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